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Abstract. - Fractional Brownian motion (fBm) has been used as a theoretical framework to 
study real time series appearing in diverse scientific fields. Because its intrinsic non-stationarity 
and long range dependence, its characterization via the Hurst parameter H requires sophisticated 
techniques that often yield ambiguous results. In this work we show that fBm series map into 
a scale free visibility graph whose degree distribution is a function of H. Concretely, it is shown 
that the exponent of the power law degree distribution depends linearly on H. This also applies 
to fractional Gaussian noises (fGn) and generic noises. Taking advantage of these facts, 
we propose a brand new methodology to quantify long range dependence in these series. Its 
reliability is confirmed with extensive numerical simulations and analytical developments. Finally, 
we illustrate this method quantifying the persistent behavior of human gait dynamics. 
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' Self-similar processes such as fractional Brownian 
ipotion (fBm) [1] arc cmTcntly used to model fractal 
phenomena of different nature, ranging from Physics or 
Piology to Economics or Engineering. To cite a few, 
fBm has been used in models of electronic delocaliza- 
tion [2], as a theoretical framework to analyze turbulence 
data [3], to describe geologic properties [4], to quantify 
qorrelations in DNA base sequences [5], to characterize 
physiological signals such as human heartbeat [6] or gait 
dynamics [7], to model economic data [8] or to describe 
network traffic [9-11]. Fractional Brownian motion Bnit) 
is a non-stationary random process with stationary 
self-similar increments (fractional Gaussian noise) that 
can be characterized by the so called Hurst exponent, 
< H < 1. The one-step memory Brownian motion is 
obtained for H = 

shows persistence and anti-persistence li H < 

While different fBm generators and estimators have been 
introduced in the last years, the community lacks consen- 
sus on which method is best suited for each case. This 
drawback comes from the fact that fBm formalism is ex- 
act in the infinite limit, i.e. when the whole infinite series 
of data is considered. However, in practice, real time se- 



whereas time series with H > ^ 



ries are finite. Accordingly, long range correlations are 
partially broken in finite series, and local dynamics cor- 
responding to a particular temporal window are overes- 
timated. The practical simulation and the estimation 
from real (finite) time series is consequently a major is- 
sue that is, hitherto, still open. An overview of different 
methodologies and comparisons can be found in [11-18] 
and references therein. Most of the preceding methods 
operate either on the time domain (e.g. Aggregate Vari- 
ance Method, Higuchi's Method, Detrended Fluctuation 
Analysis, Range Scaled Analysis, etc) or in the frequency 
or wavelet domain (Periodogram Method, Whittle Esti- 
mator, Wavelet Method). In this letter we introduce an 
alternative and radically different method, the Visibility 
Algorithm, based in graph theoretical techniques. In a 
recent paper this new tool for analyzing time series has 
been presented [19]. In short, a visibility graph is ob- 
tained from the mapping of a time series into a network 
according with the following visibility criterium: two ar- 
bitrary data (ta,ya) and (th^yh) in the time series have 
visibility, and consequently become two connected nodes 
in the associated graph, if any other data {tc, Vc) such that 
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Fig. 1: Example of a time series (20 data, depicted in the 
upper part) and the associated graph derived from the visibility 
algorithm. In the graph, every node corresponds, in the same 
order, to a series data. The visibility rays between the data 
define the links connecting nodes in the graph. 



ta <tc < h fulfills: 



yc<yb + iva - Vh) 
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In fig.l wc have represented for illustrative purposes an 
example of how a given time series maps into a visibility 
graph by means of the Visibility Algorithm. A prelim- 
inary analysis has shown that series structure is inher- 
ited in the visibility graph [19]. Accordingly, periodic se- 
ries map into regular graphs, random series into random 
graphs and fractal series into scale free graphs [20]. In 
particular, it was shown that the visibility graph obtained 
from the well-known Brownian motion has got both the 
scale-free and the small world properties [19]. Here we 
show that the visibility graphs derived from generic fBm 
series are also scale free. This robustness goes further, 
and we prove that a linear relation between the exponent 
7 of the power law degree distribution in the visibility 
graph and the Hurst exponent H of the associated fBm 
series exists. Therefore, the visibility algorithm provides 
an alternative method to compute the Hurst exponent and 
then, to characterize fBm processes. This also applies to 
fractional gaussian noise (fGn) [1] which are nothing but 
the increments of a fBm, and generic noises, enhanc- 
ing the visibility graph as a method to detect long range 
dependence in time series. 

In fig. 2 we have depicted in log-log the degree distribu- 
tion of the visibility graph associated with three artificial 
fBm series of 10^ data, namely an anti-persistent series 
with H = 0.3 (triangles), a memoryless Brownian mo- 
tion with H = 0.5 (squares) and a persistent fBm with 
H = 0.8 (circles). As can be seen, these distributions fol- 
low a power law P{k) ~ k~'^ with decreasing exponents 

70.3 > 70.5 > 70.8- 

In order to compare 7 and H appropriately, we have 
calculated the exponent of different scale free visibility 
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Fig. 2: Degree distribution of three visibility graphs, namely 
(i) triangles: extracted from a fBm series of 10^ data with 
H = 0.3, (ii) squares: extracted from a fBm series of 10"" data 
with H = 0.5, (iii) circles: extracted from a fBm series of 10^ 
data with H — 0.8. Note that distributions are not normal- 
ized. The three visibility graphs are scale- free since their degree 
distributions follow a power law P{k) ~ k ' with decreasing 
exponents 70.3 > 70.5 > 70.8- 



graphs associated with fBm artificial series of 10^ data 
with < H < 1 generated by a wavelet based algo- 
rithm [23]. Note at this point that some bias is inevitably 
present since artificial scries generators are obviously not 
exact, and consequently the nominal Hurst exponents have 
an associated error [21]. For each value of the Hurst pa- 
rameter we have thus averaged the results over 10 real- 
izations of the fBm process. We have estimated exponent 
7 in each case through Maximum Likelihood Estimation 
(MLE) [25]: 



J — 1 + n 



i=l 



(2) 



where n is total number of values taken into account, 
Xi, i = l,..,n are the measured values and Xmin corre- 
sponds to the smallest value of x for which the power law 
behavior holds. In fig. 3 we have represented the relation 
between 7 and H (black circles). As can be seen, a roughly 
linear relation holds (the dotted line represents the best 
linear fitting 7 = 3.1 — 2H). 

That fBm yields scale free visibility graphs is not that 
surprising. The most highly connected nodes (hubs) are 
the responsible for the heavy tailed degree distributions. 
Within fBm series, hubs are related to extreme values in 
the series, since a data with a very large value has typically 
a large connectivity, according to eq. 1. In order to cal- 
culate the tail of the distribution we consequently need to 
focus on the hubs, and thus calculate the probability that 
an extreme value has a degree k. Suppose that at time 
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Fig. 3: (Black dots) Numerical estimation of exponent 7 of 
the visibility graph associated with a fBm series with expo- 
nent H. In each case 7 is averaged over 10 realizations of a 
fBm series of 10'' data, in order to avoid non-stationary biases 
(the error bars are included in the dot size). The dotted line 
corresponds to the best linear fitting 7(-ff) = a — bH , where 
a = 3.1 ± 0.1 and b = 2.0 ± 0.1, and the solid line corresponds 
to the theoretical prediction "i{H) = 3 — 2H. Both results 
are consistent. Note that deviations from the theoretical law 
take place for values oi H > 0.5 and H < 0.5 (strongly corre- 
lated or anti-correlated series), where fBm generators evidence 
finite-size accuracy problems [21], these being more acute the 
more we move away from the non-correlated case H = 0.5. 

t the series reaches an extreme value (a hub) Buit) = h. 
The probability of this hub to have degree T is 

P{T) ^ Pfr{T)r{T), (3) 

where Pfr{T) provides the probability that after T time 
steps, the series returns to the same extreme value, i.e. 
B{t + T) = h (and consequently the visibility in t gets 
truncated in t + T), and r(T) is the percentage of nodes 
between t and t + T that t may see. Pfr{T) is nothing but 
the first return time distribution, which is known to scale 
as PfriT) ~ for fBm series [22]. On the other hand, 

the percentage of visible nodes between two extreme values 
is related to the roughness of the series in that basin, that 
is, to the way that a series of T time steps folds. This 
roughness is encoded in the scries standard deviation [1], 
such that intuitively, we have r{T) ~ /T = T"-^ (this 
fact has been confirmed numerically). Finally, notice that 
in this context T = fc, so eq.3 converts into 

what provides a linear relation between the exponent of 
the visibility graph degree distribution and the Hurst ex- 
ponent of the associated fBm scries: 

7(i/) = 3-2H, (5) 



Fig. 4: (Black triangles) Numerical estimation of exponent 7 
of the visibility graph associated with a noise. In each 
case 7 is averaged over 10 realizations of a series of 10^ 
data, in order to avoid non-stationary biases (the error bars are 
included in the triangle size). The straight line corresponds to 
the theoretical prediction in eq.7. 

in good agreement with our previous numerical results. 
Note in figure 3 that numerical results obtained from ar- 
tificial series deviate from the theoretical prediction for 
strongly-correlated ones {H > 0.5 ot H < 0.5). This de- 
viation is related to finite size effects in the generation of 
finite fBm series [21], and these effects arc more acute the 
more we deviate from the non-correlated case H = 0.5. 
In any case, a scatter plot of the theoretical (eq.5) versus 
the empirical estimation of j{H) provides statistical con- 
formance with a correlation coefficient c = 0.99. 
To check further the consistency of the visibility algorithm, 
an estimation of the power spectra is performed. It is well 
known that fBm has a power spectra that behaves as 1 //^, 
where the exponent (3 is related to the Hurst exponent of 
an fBm process through the well known relation [24] 

f3{H) = 1 + 2H. (6) 

Now according to eqs.5 and 6, the degree distribution of 
the visibility graph corresponding to a time series with 
noise should be again power law P{k) ~ where 

7(/3) = 4-/3. (7) 

In fig. 4 we depict (triangles) the empirical values of 7 cor- 
responding to artificial series of 10^ data with /3 rang- 
ing from 1.2 to 2.8 in steps of size 0.1 [26]. For each value 
of P we have again averaged the results over 10 realiza- 
tions and estimated P through MLE (eq.2). The straight 
line corresponds to the theoretical prediction eq.7, showing 
good agreement with the numerics. In this case, a scatter 
plot confronting theoretical versus empirical estimation of 
7(/3) also provides statistical conformance between them. 
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up to c = 0.99. 

Finally, observe that eq.6 holds for IBm processes, while 
for the increments of an fBm process, known as a frac- 
tional Gaussian noise (fGn). the relation between /? and 
H turns to be [24]: 



m) 



-1 + 2iJ, 



(8) 



where H is the Hurst exponent of the associated fBm pro- 
cess. We consequently can deduce that the relation be- 
tween 7 and H for a fGn (where fGn is a series composed 
by the increments of a fBm) is 



7(iJ) = 5-2H. 



(9) 



Notice that eq.9 can also be deduced applying the same 
heuristic arguments as for eq.5 with the change H 
H-l. 

In order to illustrate this latter case, we finally address a 
realistic and striking dynamics where long range depen- 
dence has been recently described. Gait cycle (the stride 
interval in human walking rhythm) is a physiological sig- 
nal that has been shown to display fractal dynamics and 
long range correlations in healthy young adults [27,28]. In 
the upper part of fig. 5 we have plotted to series describing 
the fluctuations of walk rhythm of a young healthy person, 
for slow pace (bottom series of 3304 points) and fast pace 
(up series of 3595 points) respectively (data available in 
www.physionet .org/physiobank/database/umwdb / [29] ) . 
In the bottom part we have represented the degree dis- 
tribution of their visibility graphs. These ones are again 
power laws with exponents 7 = 3.03 ± 0.05 for fast pace 
and 7 = 3.19 ±0.05 for slow pace (derived through MLE). 
According to eq.7, the visibility algorithm predicts that 
gait dynamics evidence behavior with (3=1 for fast 
pace, and [3 = 0.8 for slow pace, in perfect agreement with 
previous results based on a Detrended Fluctuation Anal- 
ysis [27,28]. These series record the fluctuations of walk 
rhythm (that is, the increments), so according to eq.9, the 
Hurst exponent \s H ~ \ for fast pace and 7? = 0.9 for 
slow pace, that is to say, dynamics evidences long range 
dependence (persistence) [27,28]. 
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Fig. 5: Black signal: time series of 3595 points from the stride 
interval of a healthy person in fast pace. Red signal: time 
series of 3304 points from the stride interval of a healthy person 
in slow pace. Bottom: Degree distribution of the associated 
visibility graphs (the plot is in log-log). These are power laws 
where 7 = 3.03 ± 0.05 for the fast movement (black dots) and 
7 = 3.19±0.05 for the slow movement (red dots), what provides 
P = 1 and /3 = 0.8 for fast and slow pace respectively according 
to eq.7, in agreement with previous results [27,28]. 



fGn signals (the increments of a fBm where H is defined 
as [3{H) = -1 + 2H), we have ^(7) = 
The reliability of this methodology has been confirmed 
with extensive simulations of artificial fractal series and 
real (small) series concerning gait dynamics. To our 
knowledge, this is the first method for estimation of long 
range dependence in time series based in graph theoretical 
techniques advanced so far. Some questions concerning 
its accuracy, flexibility and computational efficiency will 
be at the core of further investigations. In any case, we 
do not pretend in this work to compare its accuracy with 
other estimators, but to propose an alternative and sim- 
ple method based in completely different techniques with 
potentially broad applications. 



As a summary, the visibility graph is an algorithm that 
map a time series into a graph. In so doing, classic meth- 
ods of complex network analysis can be applied to char- 
acterize time series from a brand new viewpoint [19]. In 
this work we have pointed out how graph theory tech- 
niques can provide an alternative method to quantify long 
range dependence and fractality in time series. We have 
reported analytical and numerical evidences showing that 
the visibility graph associated to a generic fractal series 
with Hurst exponent H is, a. scale free graph, whose degree 
distribution follows a power law P{k) k^^ such that: (i) 
There is a universal relation between 7 and the exponent 
/? of its power spectrum that reads 7 = 4-/3; (ii) for fBm 
signals (where H is defined such that I3{H) = 1 + 2_ff), 
the relation between 7 and H reads i?(7) = while for 
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